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We reformulate the problem of the cancellation of the ultraviolet divergencies of the vacuum 
energy, particularly important at the cosmological level, in terms of a saturation of spectral function 
sum rules which leads to a set of conditions on the spectrum of the fundamental theory. We specialize 
the approach to both Minkowski and de Sitter space-times and investigate some examples. 
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I. INTRODUCTION 



The energy density of the vacuum, and in general 
the cosmological constant problem, has always attracted 
much attention [l], 0] . The recent discovery of cosmic ac- 
celeration [3[ has stimulated the discussion about what 
kind of dark energy could be responsible for this phe- 
nomenon [|[ ■ The cosmological constant is the most nat- 
ural candidate. The enigma arises because the zero point 
energy associated with quantum fields and naively cal- 
culated by introducing a cutoff at the Planck scale ap- 
pears to be some 118 orders of magnitude bigger than 
the observable value of the cosmological constant energy 
density. Many explanations have been put forward for 
this [5[, prominent among them are symmetry mecha- 
nisms and in particular supersymmetry. Indeed since 
the contributions to the ground state energy of Fcrmions 
and Bosons have opposite signs unbroken supersymme- 
try leads to a vanishing vacuum energy. Unfortunately 
supersymmetry, if it exists, certainly must be broken and 
at this point one appears to still be left with the cosmo- 
logical constant problem. One lesson one learns however 
is that the number of fermionic and bosonic degrees of 
freedom must be equal, if one wants to have an exact 
cancellation at least of the quartically ultraviolet diver- 
gent part of the vacuum energy 1 . 

In order to employ what one has learned from the above 
one may seek inspiration from the study of other sym- 
metry breaking, in particular chiral symmetry. Indeed if 
chiral SU (2) x SU (2) was an exact symmetry one would 
have expected the p meson to be accompanied by an ax- 
ial vector meson of the same mass, which certainly is 
not the case. A way out of this was found through the 
introduction of an asymptotic chiral symmetry [fj lead- 
ing to spectral function sum rules which are related to 
the short-distance behaviour of products of vector and 



axial- vector currents pj. The approximate saturation of 
the spectral functions with suitable combinations of par- 
ticles lead to a physically satisfactory relation between 
vector and axial- vector meson masses. 
We wish to follow an analogous procedure for the case of 
vacuum energies. Indeed we shall first obtain expressions 
for the large momentum behaviour of the vacuum energy 
in de Sitter space for arbitrary mass, spin 1/2 Fermions 
and spin and spin 1 Bosons. This is done in the next 
Section. Subsequently we shall, through the use of spec- 
tral functions, represent the contributions of an arbitrary 
number of Fermions and Bosons. On examining the con- 
tributions of all particle species to the vacuum energy 
and requiring that all divergent contributions cancel, we 
obtain constraints between diverse spectral functions and 
their moments. This of course leads to contraints on the 
mass spectra of the particles which approximately sat- 
urate the sum rules. These constraints are studied in 
Section III. Lastly our results are summarized and dis- 
cussed in the Conclusions. 



II. VACUUM ENERGY OF FREE FIELDS 

In what follows we shall consider a flat Friedmann uni- 
verse with the Robertson- Walker metric 



-dt 2 



a 2 (t)dx l dx i 



(1) 



1 The finite part of the vacuum energy (which is difficult to cal- 
culate explicitly) could be responsible for the small observable 
value of the cosmological constant 



where the cosmological radius a(t) obeys the de Sitter 
expansion law: a(t) = a^e nt . We shall calculate the 
ultraviolet divergent contributions to vacuum energy for 
free fields on this spatially flat de Sitter background. 



A. Bosonic fields 

Free bosonic field dynamics, on a Friedmann- 
Robertson- Walker (FRW) space-time, is known to be 
equivalent L 8j to that of an infinite system of time depen- 
dent harmonic oscillators (TDHO). Subsequently, at the 
quantum level, one can express the vacuum expectation 
value of the field Hamiltonian in terms of the so-called 



2 



Ermakov-Pinney variables [1, [l(J [llj] . These variables are 
introduced to define the Invariant operators one uses to 
build the Hilbert space of the solutions for the TDHO. 
Since we should be interested in the ultraviolet diver gen t 
contributions to the Bunch-Davies vacuum energy [l2j, 
we shall finally consider the modes up to an U.V. cutoff, 
and, in particular, restrict the analysis to trans-Planckian 
fc-modes of the fields, namely the limit z = Ji~ x kja ^> 1. 



1. Method of Invariants 

When the action of a quantum system can be written 
in the form of some TDHO 



S = I dtKJL- G(t) q 2 



'2 \F(t) 
with Hamiltonian given by 

H =2 P+ 2 q > 



(2) 



(3) 



and is a function of time. Setting p 2 (t ) = y/ F(t )/G(t ) 
and p(to) — 0, the invariant vacuum |0j) (which is related 
to the Bunch-Davies vacuum) coincides with the Hamil- 
tonian vacuum |0) at time to- 

Technically the full quantum evolution is determined 
once the Ermakov-Pinney equation is solved. The ex- 
act solution of the Ermakov-Pinney equation is known 
if the homogeneous, linear equation associated with it is 
solvable; otherwise the Ermakov-Pinney equation itself 
can be solved perturbatively and gives the correct adia- 
batic series. 

In particular, when handling free field evolution in dc 
Sitter space-time [13|, the adiabatic parameter should be 
chosen to be the z defined above. In the ultraviolet 
regime, one must set x = £/y/z, where £ satisfies 



a 

dz 2 



1 



H-W (z) + i 



l 



n 2 e 



(8) 



and the modified Ermakov-Pinney equation |8]) gives the 
correct adiabatic solution. 



where p — q/F, F may be interpreted as an inverse effec- 
tive mass of the harmonic oscillator while G is the mass 
multiplied by the effective frequency squared. The solu- 
tions to the Schrodinger equation are eigenstates of the 
quadratic invariant operator I defined by 



dt 



I,H 



(4) 



The quadratic invariant operator for the case ([3]) is given 
by 



/ = 



pp 



El 

F 



(5) 



where p = yF x, x is the Ermakov-Pinney variable sat- 
isfying the equation fioj ] 



fl 2 x 



(6) 



and 



— §£5-. Both the invariant §5§ 
can be written in terms of time de- 



and n 2 = FG + ^ 
the Hamiltonian (|3|) 
pendent creation-annihilation operators but, in the time 
dependent case, the Hilbert space of the solutions of the 
Schrodinger equation is generated by the linear invari- 
ant creation-annihilation operators; energy eigenstates 
are related to these solutions by a Bogoliubov rotation 
depending on time through p. F and G. The invariant 
vacuum |0j) and the invariant operator §5§ itself are not 
unique but depend on the initial conditions chosen for 
the Ermakov-Pinney equation; in any case the expecta- 
tion value of the Hamiltonian ([3]) with respect to it can 
be expressed as 



(0,|JJ|0 7 ) = i 



(7) 



2. Real scalar field 

The action for a minimally coupled, massive scalar field 
$(x, t) is 



if one considers the Fourier transform of the field 



$(x, t) = _L £ e lk ' x [fa (k, t) + *0 2 (k, t)} 
V 2v k 



(9) 



(10) 



where J dx = V and 0(k, t) is a real function, the action 
© can be rewritten as 



S = 



i=l,2 k 



dt(0i(k) 2 -w^(k) : 



(11) 



+ m 2 . We note that one can elimi- 



where lu 2 

nate the overall factor 1/2 in (|TT|) by observing that 
0i (k) = (j>i(—k) and d> 2 (\s.) = — 2 (— k) thus restrict- 
ing the sum over the subspace k+ of the independent 
degrees of freedom. For each independent mode k, i, one 
then obtains the following Hamiltonian 



2a 3 2 



(12) 



where 7Tj(k) = a 3 0i(k); the full Hamiltonian can then be 
written as Hs = J2i=i i Sk+ ^O*)- Finally, to calculate 
the vacuum expectation value ^ , in the continuum limit 



v 4^ 



(2tt) 3 



k 2 dk dCl, 



(13) 



3 



it is sufficient to just keep 
complete solid angle 



bi and integrate over the 



(0i\H s \0i) V- 



V 

1 

2^ 



k 2 (0i 



^ J k 2 E s Q {k)dk 



Tr(k) 2 a^im 

2a 3 2 



|0/)dfc. (14) 



For the case of a scalar field, one has, for a mode k, 
Fg = a~ 3 , G$ — a3uj k and the Ermakov-Pinney equation 



can be solved exactly with 



xs 



1 nH^^H^jz) 
2H 



(15) 



where the Hu (z) are the Hankel functions and v — 
! — Sr< We note that the solution (fT5|) has been 
chosen to fulfil the requirement that |0/) coincides with 
Bunch-Davies vacuum. One can finally calculate Eo(k) 
exactly and evaluate it in the z » 1 limit 



Z a 



4fc 



m 2 (2H 2 - rn 2 ) a 3 f 1 
16P + ° 



(16) 



3. Real Vector Field 

The Proca action for a minimally coupled, massive vec- 
tor field (x, t) = (Aq , A) is given by 



S\A b 



(Fx, 



1 



F^F^ 



M 2 



A»A„ 



(17) 



where F^ = V^A V ~ V„j4 m = d M A„ - (with V M 

the covariant derivative). We note that in the massive 
case the vector field satisfies the Lorentz gauge condition 
VM M = and has three independent physical degrees 
of freedom; in the massless case, instead, to the Lorentz 
condition one should add Aq = 0, leading to V ■ A = 
(radiation gauge) and A^ reduces to two independent 
physical components. The final expression one obtains 
for the action (|17p . without any gauge fixing, is 



S[A, 



VAA 



.4 



V I,] 2 - 2A- (VA 



M 2 (a 2 Al - A 2 ) 



(18) 



For the massive case, the variation of the action with 
respect to Aq leads to the Gauss constraint: 



1 



V ■ A - V Ai 



M 2 A Q = 0. 



(19) 



If one considers the Fourier transform of the field com- 
ponents 



ikx 



(20) 



then the action can be expressed in terms of a^(k) 
(ao(k), a(k)) as 

+ a 2 w 2 a (k)a (-k) 



k A a(k) • k A a(-k) 2 2 



M 2 a(k) • a(-k) 



(21) 



with w 2 = M 2 + ^j- and the Gauss equation IT9|) for each 
mode is 



a (k) 



a(k). 



(22) 



The action (|21|) can now be written in terms of the com- 
ponents <3(k). If one decomposes <3(k) into the sum of a 
longitudinal vector 

d L (k) = (d?(k) • is) k = [«l (k) + i u £ (k)] fc (23) 
V2 



with k = k/vk 2 , plus a transverse vector 
a T (k) = a(k) - a L (k) 

= 4 S [t4°(k) + i^(k)]gi(k) (24) 

V ^ £=1,2 

where the e^(k) form an orthonormal basis, e^(k)-ej(k) = 
5ij, for the transverse space, ej(k) ■ k = 0, the action (|2~1]) 
can finally be written as 



k+ 



S L [u L (k)]+S L [v L (k)} 
]T (s T [u$(k)} + S T [vP(k)] 



i=l,2 



where 



S L [p(k)] = M 2 | dt| {^m 2 ~ ^(k) ; 



and 



5 T [^(k)]= / dt-^k) 2 -^^) 2 ) 



(25) 



(26) 



(27) 



We note that, owing to the overall M 2 factor in front 
of (|26p . one recovers the correct physical limit for the 
massless case: in fact only the transverse contributions 
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remain when M — > 0. In the latter case, the full Hamilto- 
nian vacuum expectation value (v.e.v.), in the continuum 
limit (fT3|) . is given by 



(0l\H T \0l) V^oo 1 



V 



k 2 E> {k)dk 



k 2 (0j 



7r T (k) 2 aulMV 2 



2a 



|0j)dfe (28) 



where 7r:r(k) = a^r(k), the factor 2 accounts for the con- 
tributions coming from different transverse polarizations, 
and the M — > limit should be taken. In the massive 
case one has also the longitudinal part 



(0i\H L \0i 



V 

1 

2^2 



K->oo 1 

* 2^2 



k 2 E^(k)dk 



k 2 {Qi\ 



^(k) 2 , oM^ £ (k) s 
2aM 2 



+ 



|0/)dfc(29) 



B. Fermionic Field 

The action for a massive spin \ field ^(x, i) is: 



-At*^l (34) 



where V M is the gauge covariant derivative for spin 1/2 
fields, 7 M = e a M (x)7 a and eJ J '(x) are the vierbein. In 
FRW space-time, on Fourier expanding 



(35) 



the action can be rewritten in a compact formfl^ 



as 



■ dt - M(k) 



V-(k) (36) 



With 7T L (k)= -^(k). 

The Ermakov-Pinney equation ^ can be solved exactly 
for the transverse part since Ft = a -1 , Gt — auj 2 , and 



2H 



(30) 



with v= J\- jjr is the solution which is associated 
with the Bunch-Davies vacuum. The vacuum energy, in 
the large z limit, is given by 

T z »i k , M 2 a M 4 a 3 
E (k) 



2 a 



4fc 



16/c 3 



1 



(31) 



We note that quadratic and logaritmic divergent terms 
disappear in the massless case. 

The longitudinal field equations are obtained by setting 



F L = 



aM 



-, Gl = aM 2 and cannot be solved exactly; 



however one can estimate the vacuum energy of (|29[) in 
the adiabatic, z >> 1 limit by expanding the solution of 
Eq. (J8]) as a power series in z 

z »i 1 I". . 2H 2 -M 2 



XL 



1 

5 (M 



1 



m 2 z 2 

12M 2 H 2 



4H 4 ) 



32H 4 z 4 

and then expanding ([7]) to obtain 



(32) 



2>1 



(H 2 + M 2 



2a 4fc 
M 2 (6H 2 + M 2 )a 3 
16~P 



3 . (33) 



Interestingly enough, the longitudinal vacuum energy is 
slightlty different from that of a scalar field because of 
the curved space kinematics. 



where M(k) is the 4x4 matrix 



M(k) 



u a ■ — 

' a 



(37) 



and a are the standard 2x2 Pauli matrices. The Euler- 
Lagrange equation obtained for a spinor ^(k) is 



iip(k) = M(k)V>(k) 



(38) 



and the Hamiltonian for the system can be easily calcu- 
lated in terms of ip(k) as 

if / = ^if / (k) = ^V(k) t M(k)V(k). (39) 



The general solution of the equation of motion (|3"8"|) can 
be written as a superposition of the 4 independent solu- 
tions 



with r = 1, 2 and 



iHl l \z) 



%H$[z) 



(40) 



(41) 



with r = 3,4, where 2/ = ^ — i£, Nf — exp [7r fJ./'H], 
r = 1, 3(2, 4) correspond to spin up (down) for the Pauli 

spinors an d the depend on time through z. 

Note that the independent solutions (I40l4ip reduce to the 
usual static solutions in the z — > +oo limit and they may 
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then be associated with the time-independent creation- 
annihilation operators b r (k), b r (ky, d r (k), d r (ky . One 
defines the Bunch-Davies vacuum \0bd) by 



b r (k)\0 BD ) = d r (k)\0 BD ) =0 



(42) 



with 



{S r (k),S r (k)t} = (dr(k),d r (k^\ = 1. (43) 



In terms of the above operators the Heisenberg field op- 
erator ip(k, t) can be written as 



r=l,2 

+(-ir +i 4 r) (k,i)d r (-k)t 



(44) 



and the vacuum expectation value of the Hamiltonian 
([55]) can be evaluated in the continuum limit (fTS")) 



1 



r 



2vr 2 



k 2 dk (0 BD \H f (k)\0 BD ) 



(45) 



the Einstein equations coming from different free fields 
using the Kallen-Lehmann spectral functions formalism. 
Some specific models will then be analyzed in detail. For 
the Minkowski case one can obtain some general results 
which will later be generalized in a non trivial way to the 
de Sitter case. 

A. Spectral function general formalism 



Starting from a semiclassical approach we consider the 
Friedmann equation for the homogeneous scale factor 
a(t) in J!]) 



3^2 _ 87rG {^BD\Hj\Q B D 



- 

3 



V 



(48) 



where is the sum over all bosonic and fermionic field 
(representing the particle content evolving on the space- 
time manifold) Hamiltonians averaged with respect to 
the corresponding Bunch-Davies vacua. The spectral 
function pa(x) can be introduced if we replace the above 
sum by 



by using (I42ti44|) . We note that each i?/(k) contains the 
contributions both of a particle and an antiparticle with 
opposite spin and momentum; in terms of the Hankel 
functions, one obtains the following exact expression: 

(0 BD \H f (k)\0 BD ) = 2 [m (h&H™-! ~ HPH$> 



and, on counting a single particle/antiparticle degree of 
freedom, in the large z limit, one finds 



i Q 



2a 4k 



16 fc 3 



(47) 



III. SUM RULES AND CANCELLATION OF 
ULTRAVIOLET DIVERGENCES 

In this section we consider the general condition for 
the cancellation of divergent vacuum contributions to 



E 



2tt 2 



k 2 dk 



dm 2 ps{m 2 ) 



+ J dM 2 (p v (M 2 ) - v S(M 2 )) 
+ J dM 2 2p v (M 2 ) 
+ [ dp 2 2p F (p 2 ) 



(49) 



where Pa(%) have support for < x < oo, vq8{M 2 ) is the 
number density of massless vector particles, and Dirac 
spinors account for double the contribution of Majorana 
and Weyl particles. Moreover the second line of Eq. (|4"9"|) 
counts the longitudinal vector degrees of freedom while 
the third line gives the transverse contributions. On a 
de Sitter background, setting m 2 = x with A being the 
ultraviolet cutoff (k < aHA) we adopt to regularize the 
divergent integral L dk , Eq. (|48]) can be finally rewrit- 
ten as 



6 



^ = H I z 2 dz I dx[ps{x)E^{x,z) + {p v (x)-v 6{x))E^{x 1 z) + 2p v (x)E^(x,z) + 2p F {x)E^{x,z)] 



dx[A 4 ^-Fi(x) + ^F 2 (x) 



8 



(50) 



where 

Fi(x) = p s (x) + 3p v (x) - v S(x) - 2p F (x) , (51) 
F»(x)= {U 2 + x) p s (a?) + (3 x + H 2 ) pv (a?) 

-v (x + H 2 )5(x)-2xp F {x), (52) 
F 3 {x) = x (2H 2 - x) ps(x) - x (6H 2 + 3x) p v (x) 

+v Q x (6H 2 + x) 5{x) + 2 x (H 2 + x) p F {x). (53) 

In order to obtain a finite contribution on the right hand 
side of Eq. 1[5D]). the spectral functions describing the 
particle content of the universe should be such as to make 
(15H53[) vanish simultaneously and this occurs when 



J dxFi{x) = 



(54) 



dxxFi(x) — Ti / dx \vq8(x) — ps{x) — pv{x) \ (55) 



dxx Fi{x) — 2H / dx x [ ps(x) — 3 pv(x) + 3vq5(x) 



+p F (x) 



(56) 



In the H — > limit the conditions (|54ti56[) reduce to the 
very compact form 



dxx l F 1 (x) = 0, i = 0,1,2. 



(57) 



We note that the expression (|50j) is completely general 
and refers to arbitrary mass distributions. In the follow- 
ing sections we shall analyze models wherein the mass 
distributions are discrete and the Pi{x) are just super- 
positions of Dirac delta functions multiplied by integer 
coefficients, pa — *52i n i{A)8(x — Xi), and describe di- 
verse particle multiplets. 

For % = Eq. (fBT)) gives the condition of the equality 
of numbers of bosonic and fermionic degrees of freedom. 
For i = 1, 2 one has 



which can be regarded as hyperplane and hypersphcrc 
equations respectively in term of squared mass variables, 
provided the bosonic or fermionic field content is fixed. 



B. Hyper sphere and hyperplane 

While the cancellation of the quartic divergences re- 
quires the equality of the numbers of fermionic and 
bosonic degrees of freedom, the conditions for the cancel- 
lation of quadratic and logarithmic divergences are more 
involved. However, for the Minkowski case these condi- 
tions can be represented as hyperplane and hypersphere 
equations respectively. 

Let us consider the hypersphere of dimensionality n—1 
embedded into the space of n dimensions, given by the 
equation 



and the hyperplane given by the equation 



Vi = A ' 



(58) 



(59) 



We look for the points of intersection between the hy- 
persphere and hyperplane such that all the coordinates 
Hi have non-negative values. Then the minimum value 
of the parameter A is A = 1 which corresponds to the 
case when only one of coordinates is non-zero, while the 
maximum value of A is A = y/n, when all the coordinates 
are equal (yi = for Vi). Thus 



1 < A < y/n. 



(60) 



The intersection between the hypersphere and the 
hyperplane f|59[) is a (n — 2)-dimensional hypersphere, 
whose center has equal coordinates yi for all i and ra- 
dius yl — X 2 /n . All the points of the hypersphere of 
intersection 5^ 1-2 ' can be represented by the vector 



1 A2 - 
1 u, 

n 



where 



€T= (!,-•• ,1), 



u ■ u = 1, u ■ v = 0. 



(61) 

(62) 
(63) 
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The coordinates of the points, belonging to S^ n 2 - ) are 

(64) 



A / A 2 

Vi = - + V 1 Ui 

n V n 



Let us find the minimum and maximum possible values 
of Ui. Without losing generality we will focus on j = 1. 
We shall look for a solution in the form of a vector 



= {<*,&), 



(65) 



where f3 is a (n — l)-dimensional vector. The conditions 
([63]) give 



and 



a 2 + f] ■ (3 = 1 



a + P ■ I = 0, 



(66) 



(67) 



where I is a (n — l)-dimensional vector with all the co- 
ordinates equal to 1. Then 



0-0 ((n-l)cos 2 x + l) = 1. 



(68) 



where \ is an angle between the vectors f3 and /. Appar- 
ently, a 2 has maximum value when the angle cosx = ±1. 
Then 



a = ± 



n- 1 



(69) 



which is non-negative for A > 1 and hence is always non- 
negative for our range of A. 

Let us now consider the case A < \Jn — 1 . In this 
case some of the points of g^ 1-2 ) have negative values 
of coordinates and should be excluded from considera- 
tion. These patches of the hypersphere have the form of 
n spherical caps, whose pole points are the points where 
one of n coordinates has its minimum and, hence, for 
A < \/n — 1 negative value. The angle characterizing the 
size of these spherical caps (i.e. the angle between the 
unit vector pointing to the pole and the vector pointing 
to the intesection of S^™ -2 ) with the hyperplane yi = 0) 
is 



A 



= arccos ■ 



V(n-A 2 )(n-1) 



(75) 



One can find also the angle between two vectors pointing 
to two different "poles" : 



X : 



arccos — 



1 



1 



(76) 



In the case for which 29 > \ two spherical caps can 
intersect each other. A simple calculation shows that it 
occurs if A < y/n — 2. Let us note that this does not 
mean that acceptable solutions do not exist. Indeed, as 
we explained above, when one of the coordinates has 
the maximum value (|7ip all the other coordinates have 
non-negative values ([74")) : for the case A > 1 these values 
are positive. 



(3 = T 



y/n(n- 1) 



(70) 



Substituting the maximum and minimum values of a 
from Eq. (|69[) into Eq. (|64p one can find the maximum 
and minimum values of one of the coordinates yf. 



A 



1 - 



A 2 n-1 



Vi 



(71) 



(72) 



We note that when one of the coordinates takes an ex- 
ternum value, all the other coordinates have equal values 
(see Eq. ([70]) ). Now, on using Eq. l[72| we can get the 
condition for the positivity of all the coordinates yi of 
the points lying on the (n — 2) - dimensional sphere. Re- 
quiring yi m i n > we have 



A > \Jn - 1. 



(73) 



Curiously, when one of coordinates acquires the maxi- 
mum value (|71|) all the other coordinates have value 



A / A 2 
Vi = Z-\ l - 



1 



n y/n(n - 1) ' 



(74) 



C. Some simple models 

In our case the parameter A is nothing more than 



A 



(77) 



where rriio and to^m are the masses of Dirac and Majo- 
rana spinors respectively, while D and M are the num- 
bers of these spinors. The number n is the total number 
of fermionic degrees of freedom 



n = AD + 2M + 2W, 



(78) 



where W is the number of Weyl spinors. The number of 
bosonic degrees of freedom should be equal to the number 
of fermionic degrees of freedom n and the boson masses 



(squared) are given by yt J 4 J2iLi TO ? 



It 



"iD ~ " "HM 

is easy to see that the value of A cannot take the minimum 
value A = 1 (cf. Eq. ([60])). 

Let us consider the model with two degrees of freedom 
(n = 2). In this case we have a Majorana spinor and the 
value of the parameter A is A = \[2 and the circumfer- 
ence and straight line have the only point of intersection, 
when yi = j/2- Obviously this situation corresponds to 
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to two scalars with the same mass, which coincides with 
that of the Majorana spinor. Therefore in this case the 
conditions to have the required cancellations are so se- 
vere that the only possible solution is the one which can 
be obtained by imposing supersymmetry for bosons and 
fermions. We shall see that for larger systems these con- 
straints are not so stringent. If instead of Majorana one 
has a Weyl spinor, the bosonic part of the spectrum will 
be represented by two massless scalar fields or by one 
electromagnetic field. 

The model with the four degrees of freedom repre- 
sented by one Dirac spinor or by two Majorana spinors 
with identical masses (n = 4, A = 2) also has only a triv- 
ial solution: four scalar fields with the same mass or one 
scalar field and one massive vector field whose masses 
coincide. The case with two Majorana spinors with dif- 
ferent masses is more interesting. In this case n = 4, 
and 



case with one scalar and one vector massive fields cor- 
responds to the vector u, considered in the preceding 
section, having an extremum value component while all 
the others have equal values. Thus, when the condition 
([73)) is satisfied there are two solutions for the system 
of equations (|8ip and (|82p. corresponding to the scalar 
field having maximum and minimum values, while if this 
condition is not valid only one solution survives. 

Let us now consider the model with 6 degrees of free- 
dom, wherein the fermion sector is represented by one 
Dirac and one Majorana spinors. To begin with, con- 
sider the bosonic sector containing two massive vector 
fields. The sum rules in this case lead to 



3a;i + 3x 2 = 4yi + 2y 2 , 



3^ + 3^ = 4^ + 2^. 



(85) 



(86) 



A = V2 



\Jm\ + ' 



(79) 



Now on using the results of the preceding section it 
is easy to see that the prohibited caps do not exist if 
A > \/3 (see Eq. (|T3"|) ). Substituting this condition into 
Eq. (fT9"]l one re-expresses this condition in terms of the 
relation between two fermion masses: 



2 + %/3 > mf/m 2 > 2 — v/3. 



(80) 



If this condition is violated the prohibited caps exist, but 
they do not intersect each other. 

In the following we shall investigate some specific mod- 
els by directly studying the conditions (|55| and (|56|) and 
using a slightly different notation: xi and yi will now 
denote the square of the boson and fermion masses re- 
spectively. 

Let us consider a particular case with two fields in the 
bosonic sector: one scalar and one vector massive fields 
are present. Now, the sum rules lead to 



where y\ and y 2 represent the Dirac and Majorana 
fermion masses squared. The solution of equations (|85|) 
and ([Ml) is 



£1,2 



2yi + 2/2 ± - Vi\ 



(87) 



This solution is positive for both x\ and X2 provided the 
following condition is satisfied: 



— < 4 + 3^2. 
2/i 



(88) 



Let us now consider the bosonic sector which contains a 
massive vector field, an clectromagentic field and a scalar 
field. In this case the sum rules are 



3.ti + x 2 = 4yi + 2y 2 , 



3xl + x 2 = Ayl + 2y 2 



(89) 



(90) 



3xi + x 2 = 2j/i + 2y 2 , 



(81) The existence of positive solutions for x\ and x 2 depends 
on the relation between fermion masses. Namely, if 



3xi 



4 = 2yf + 2y 2 2l 



(82) 



where x\ and x 2 are the vector boson and scalar masses 
squared, while y\ and y 2 are the Majorana spinor masses 
squared (y 2 >yi)- On solving Eqs. (|8l"j) and (|82j) we get 



(83) 



2/i+ J/2 , V3 
x\ = — ^ — ± ~~Q~( y2 ~ Vl > 



yi < 3V2-4 
2/2 - 2 



there is one solution 



.2/2 y 2 2/12/2 

xi=yi+ Y-y-2-— ■■ 



(91) 



(92) 



x 2 



2/1 +2/2 



V3, 

+ —\v? 



2/i ) 



(84) 



Let us note that the first solution for x 2 is negative if 
Hi- = < 2 — \/3- Obviously this condition coincides 
with the condition ([SO]) . This is quite natural, since the 



x 2 = 2/1 



If 



2/2 



yj 

12 



2/12/2 



3V2-4 < m < l 
2 - y 2 - 4 



(93) 



(94) 
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then two solutions exist. One of them is the solution 
([HI]), ([Ml) and the second solution is 



xi = Vx + — + 



yj 

12 



2/12/2 



.2/2 „ Ivl 2/12/2 

• T2 = yi + y- 3 Vi2-— - 



(95) 



(96) 



Finally, if — > | solutions do not exist. 



D. Simple models in the presence of the 
cosmological constant 

For the case of the cosmological constant being differ- 
ent from zero, the relations between masses of the fields, 
leading to the cancellation of divergences in the vacuum 
energy become more involved. 

Consider the case of two degrees of freedom. In this 
case the fermionic sector is represented by one Majo- 
rana spinor with mass squared denoted by y, while the 
bosonic sector is represented by two scalar fields with 
masses squared denoted by x\ and x 2 The sum rules are 



x x + x 2 = 2y - 2H 2 



2H 2 (x x 



- x 2 ) 



2y(y + H 2 



(97) 



(98) 



Non-negative solutions of this system of equations are 
possible provided the Majorana fermion mass satisfies 
the constraint 



y> 



7 + V33 



(99) 



In this case the masses of the scalar fields are 

xx,2 = y-H 2 ± \JhyH 2 - 3H 4 . (100) 

For the case of 4 degrees of freedom, one can have as 
a non-trivial solution for the fermionic sector a Dirac 
fermion. If the bosonic sector is represented by an elec- 
tromagnetic field and two scalar fields, the sum rules are: 



xx + x 2 = 4y - 2H 2 



x\ + xl- 2H 2 {xx + x 2 ) = Ay{y + H 2 ) . 



(101) 



(102) 



The system of equations (|101[) and (|102[) has non trivial 
non-negative solutions and the expected condition x\ = 
x 2 = y = is recovered as Tl tends to zero. In the 
presence of non zero 7i if 



2H 2 <y<-[,> 



19 )n 2 



one can easily obtain 



351,2 = 2y - H 2 ± y/lOyH 2 - 2y 2 - 3H 4 



(103) 



(104) 



where the lower bound in (|103|) is required for x\ t2 to be 
positive while the upper bound is needed for the reality 
of the solutions. 

The last case which we treat here is the bosonic sector 
represented by a scalar and a vector massive fields. The 
sum rules take the following form; 

3si +x 2 + 2H 2 = 4y, (105) 



3x 2 + x 2 2 + 2H 2 (3xx - x 2 ) = Ay{y + H 2 ). (106) 

If the condition y > 2TL 2 is satisfied there are two non- 
negative solutions: 



and 



xx 



x 2 



Xx 



y-H 2 



H 2 (H 2 +y) 



y + H 2 -J H2{H2+y) 



n2 _.lH 2 (H 2 +y) 



x 2 = y 



n 2 



H 2 (H 2 +y) 



(107) 



(108) 



(109) 



(110) 



In the limiting case y = 2Tt 2 only the solution (|107p 
and p08p is valid because in formula (|109|) the mass of 
the massive vector boson vanishes and this breaks the 
equality of the boson and fermion degrees of freedom. 

As a consequence of our approach it appears that the 
particle mass spectrum and the Hubble constant of the 
De Sitter space are related. Of course one may attempt 
to introduce some dynamics and an evolving mass spec- 
trum, which however we expect to still be connected to 
the Hubble parameter. We plan to return to this in the 
future. 



IV. CONCLUSIONS 

In this paper we considered the problem of vacuum 
energy in quantum field theory and cosmology. This is 
often associated with the so called cosmological constant 
problem. Indeed, the contribution of quantum vacuum 
fluctuations to the energy-momentum tensor behaves as 
a cosmological constant [l|. However, the cosmological 
constant, in principle, can also be of non-quantum ori- 
gin and be present in the theory as one of fundamen- 
tal constants. The real vacuum energy problem consists 
in the absence of a well established and justified pro- 
cedure for the renormalization of ultraviolet divergences 
in the energy-momentum tensor analogous to the the- 
ory of renormalization in standard quantum field the- 
ory in a flat (absence of gravity) space-time. Indeed, in 
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quantum field theory the ultraviolet divergences are ab- 
sorbed in the renormalization of a finite number of mea- 
surable constants, or, in other words, all the observable 
quantities become finite due to the introduction of in- 
finite counterterms into the bare Langrangian (see e.g. 
[l5L Il6l | ) . As far the very strong ultraviolet divergences 
in the energy-momentum tensor are concerned they are 
set equal to zero by the normal or Wick quantization 
of quantum fields. The last step could be justified by 
the fact that one always measures the differences be- 
tween energy levels and not the absolute values of energy. 
However, this justification fails in the presence of grav- 
ity, because of the very structure of the Einstein equa- 
tions, which connects the curvature of spacetime to the 
energy-momentum tensor. Thus, one should also take 
into account the contribution of the vacuum expectation 
value of the energy-momentum tensor on the right-hand 
side of Einstein equations. Indeed one may study mod- 
els wherein vacuum oscillations drive some stage of the 
expansion of patches of the universe can be studied (see 

El). 

It is well known that the naive calculation of the vac- 
uum energy of quantum fluctuations using a cutoff on 
the Planck scale gives a huge value for it. One has 
two possible ways out of this situation. The first one is 
the construction of a consistent renormalization theory 
for the ultraviolet divergences of the energy-momentum 
tensor, which is a very difficult task, since we do not 
have such relatively simple tools as the renormalization 
of some known physical quantities as, for example, in 
quantum electrodynamics (see, however, 18], where an 
attempt to fix the renormalization of an effective cos- 
mological constant based on some self-consistency con- 
ditions was undertaken). On the other hand one may 
require the exact cancellation of ultraviolet divergences 
in the energy-momentum tensor. This is the approach 
attempted in the present paper (a good review of anal- 
ogous approaches is given in [231 ] . further we also wish 
to mention a related flat space approach in the context 
of induced gravity [241]). Our idea can be simply formu- 
lated as follows: if you do not know what to do with 
these ultraviolet divergences, just try to eliminate them 
by introducing the condition of their cancellation as a 
quantum consistency condition of the theory. Indeed, the 
cancellation conditions on the spectral functions man- 
tioned in the Introduction and studied in the third sec- 
tion of this paper is one of the oldest examples of such 
approach. Other well-known examples of quantum self- 
consistency are connected with conditions for the cancel- 
lation of quantum anomalies. A classical example is the 
mechanism suggested by Glashow, Illioloulos and Maiani 
[l!|, who introduced the fourth quark into the theory of 
electroweak interactions to suppress the chiral anomaly. 
Later this quark was discovered experimentally and was 
called "charmed" . Another spectacular example is the 
appearence of critical dimensions in string and super- 
string theories [20]. The theory of a bosonic string can 
be formulated consistently in a spacetime of 26 dimen- 



sions, while the critical dimensionality for superstrings is 
equal to 10. Analogous relations between the dimension- 
ality of spacetime and its matter content arise also in the 
application of the BFV-BRST quantization mechanism 
to quantum cosmology [2 If. Yet another restriction on 
the particle spectrum of a theory based on the require- 
ment of the normalizability of the wave function of the 
universe was obtained in [22| . 

Thus, our work, which uses the condition of cancella- 
tion of ultraviolet divergences in the vacuum energy to 
arrive to some special conditions on the spectrum of the 
fundamental theory, is in kinship with a rather fruitful 
and ramified stream of the development of modern quan- 
tum theory. Let us now mention what are the similar- 
ities and the differences between our approach and the 
supersymmctry one. In the case of exactly supersym- 
metric models in flat spacetime there is an exact cancel- 
lation of the vacuum energy and not only of its divergent 
part. This cancellation occurs because the fermion and 
boson contributions to the vacuum energy enter with op- 
posite signs. However, exact supersymmetry cannot be 
implemented for the construction of a realistic theory of 
elementary particles. Moreover there are serious diffi- 
culties arising in the formulation of supersymmetry in 
curved background spacetimes. Here we share one com- 
mon feature with supersymmetry: the equal number of 
fermionic and bosonic degrees of freedom, which is in- 
dispensable for the cancellation of quartic ultraviolet di- 
vergences. The general conditions for the cancellation 
of quadratic and logarithmic divergences are much more 
flexible than the requiring of exact supersymmetry. On 
analysing some simple examples we have seen that there 
are many opportunities of satisfying these conditions. At 
the same time we do not require the cancellation of the fi- 
nite part of the vacuum energy. It could be different from 
zero and be responsible for the observable cosmological 
constant. 

We note that the sum rule constraints may generally 
allow for theories whose particle spectrum fit the con- 
tent of the stardard model in the low mass sector. It is 
widely believed that some kind of extension of the stan- 
dard model is needed and many of them, such as the su- 
persymmetric ones, have indeed been considered. There- 
fore our approach may be helpful to encode the minimum 
possible number of constraints on the spectrum of the 
theory necessary for the cancellation of the divergences 
of the vacuum energy, even in the presence of an effec- 
tive cosmological constant. So, one may think of some 
kind of grand unification theory, where sum rules can 
be realized. Perhaps, such a theory should contain not 
only additional supermassive gauge bosons but also some 
additional families of fcrmions. A straightforward exten- 
sion of our analysis can be done for the case of extra 
spatial dimensions, one has simply to take into account 
the dimension dependent number of physical boson and 
fermion degrees of freedom and write the general identity 
necessary to saturate the corresponding spectral function 
sum rules. We believe that all these topics deserve further 
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investigation, not to mention the introduction of higher We wish to thank G. E. Volovik for useful correspon- 

spins. dence. 
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